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Earlier  papers  by  Murty  [16]  and  Fathi  [7]  have  exhibited 
classes  of  linear  complementarity  problems  for  which  the  computational 
effort  (number  of  pivot  steps)  required  to  solve  them  by  Lemke's 
algorithm  [13]  or  Murty' s algorithm  [15]  grows  exponentially  with  the 
problem  size  (number  of  variables).  In  this  paper  we  consider  the 
sequences  of  complementary  bases  that  arise  as  these  problems  are 
solved  by  the  aforementioned  algorithms.  There  is  a natural  correspond- 
ence between  these  bases  and  binary  n-vectors  through  which  the  basis 
sequences  can  be  identified  with  particular  hamiltonian  paths  on  the 
unit  n-cube  and  with  the  binary  Gray  code  representations  of  the 
integers  from  0 to  2n  - 1. 
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OBSERVATIONS  ON  A CLASS  OF  NASTY 
LINEAR  COMPLEMENTARITY  PROBLEMS 

by 

Richard  W.  Cottle 

1 . Introduction. 

From  the  theoretical  standpoint,  the  "complementary  pivot  methods" 
of  Lemke  [13],  Murty  [15],  and  Dantzig-Cottle  [5],  will  always  solve 
the  nondegenerate  linear  complementarity  problem  (q,  M)  of  order  n : 
find  w and  z satisfying 


Iw  - Mz  = q 
w > 0 , z > 0 
zTw  = 0 

when  M e Rn*n  is  a P-matrix  and  q s Rn  is  arbitrary.  Empirically 
it  has  been  found  that  these  methods  work  quite  efficiently,  often  solv- 
ing the  problem  by  roughly  0(n)  pivots.  Encouraging  though  this  may 
be,  it  need  not  always  be  the  case.  Indeed,  Murty  [16]  and  Fathi  [7] 
have  each  exhibited  classes  of  linear  complementarity  problems  (indexed 
by  the  order,  n ) requiring  either  2n  or  2n  - 1 pivots  depending 
on  which  of  two  complementary  pivot  methods  is  used  to  obtain  the  solution. 

These  contributions  to  the  subject  of  computational  complexity  of 
algorithms  deal  with  specific  numerical  problems  (q,  M)  . In  Murty's 
paper,  M is  a lower  triangular  P-matrix  with  positive  entries  below 
the  main  diagonal,  and  hence  it  is  necessarily  nonsymmetric.  In  Fathi's 
paper,  M is  symmetric,  positive  definite  and  positive.  As  such,  these 


matrices  (to  be  specified  more  precisely  in  Section  4)  look  rather  nice 
insofar  as  the  constructive  existence  theory  is  concerned. 


The  primary  purpose  of  this  paper  is  to  point  out  that  there  is 
something  quite  special  about  the  sequences  of  bases  through  which  the 
algorithms  run  when  applied  to  the  numerical  examples  of  Murty  and  Fathi. 
A secondary  purpose  of  the  paper  is  to  characterize  the  subclass  of  P- 
matrices  and  vectors  q for  which  a particular  pivoting  sequence  is 
followed.  The  necessary  and  sufficient  conditions  given  here  are  ex- 
pressed in  terms  of  determinantal  inequalities;  other  equivalent  formu- 
lations are  possible. 

The  plan  of  the  paper  is  the  following.  In  Section  2,  we  intro- 
duce some  notation  and  background  especially  relevant  to  the  main  point. 
In  Section  3,  we  explore  some  properties  of  a particular  function  related 
to  the  special  pivot  sequence  arising  in  the  Murty-Fathi  examples.  In 
Section  4 we  analyze  the  aforementioned  examples  and  show  that  they  do 
indeed  give  rise  to  the  pivot  sequence  we  have  in  mind.  In  Section  5, 
we  establish  the  relationship  between  the  data  and  the  pivot  sequence 
followed  by  the  complementary  pivot  methods  of  Lemke  and  Murty.  Finally, 
in  Section  6,  we  mention  some  interesting  related  literature. 
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2.  Notation  and  background. 


For  the  sake  of  brevity,  we  omit  most  of  the  customary  preliminary 
material  and  refer  the  reader  to  Murty  [16]  and  Fathi  [7]  for  the  stan- 
dard concepts  and  notations  of  linear  complementarity  theory.  To  some 
extent,  we  rely  on  these  papers  for  accounts  of  the  "complementary  pivot 
methods"  to  be  discussed  here. 

However,  there  is  one  item  so  essential  to  our  development  that  it 
must  be  reviewed.  Obviously,  the  equation 

Iw  - Mz  = q (1) 

is  central  to  the  study  of  (q,  M)  . One  is  interested  in  complementary 
bases  in  the  matrix  [I,  -M]  . These  are  nonsingular  matrices 
B = [B  -j , ....  B>p]  such  that 

B . e {I  -M  •}  , j = 1 , ....  n , (2) 

•J  *J  *J 

i.e.,  the  j-th  column  of  B is  either  the  j-th  column  of  I or  the 
j-th  column  of  -M  . Such  matrices  give  rise  to  complementary  basic 
solutions  of  (1).  A nonneqative  complementary  basic  solution  of  (1) 
solves  (q,  M)  . 

Given  a complementary  basis  B , the  set 

pos  B = {q:  q = Bv  , v _>  0} 

is  called  the  complementary  cone  relative  to  B . 

Again,  suppose  a complementary  basis  B is  given.  To  this  basis 
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we  associate  a binary  row  vector  c e R as  follows.  For 
j = 1 , . . . , n let 


c . = 
J 


0 if  B = I . 

•J  *J 

1 if  B . = -M  . . 

•J  *J 


(3) 


When  M e P (the  class  of  matrices  with  positive  principal  minors), 
we  must  have  I . t -M  . , so  there  is  absolutely  no  ambiguity  in  de- 

J J 

fining  c.  as  in  (3).  In  other  cases,  one  could  use  the  definition 

J 


Cj 


0 if  B..  = [I,  -M].  . 

1 if  B.j  = Cl-  -«.n+J  ‘ 


(3') 


According  to  (3),  the  vector  corresponding  to  the  identity  matrix  I 
is  (0,  ....  0)  e Rn  . 

As  is  well  known,  there  are  2n  binary  vectors  in  Rn  , and  when 
Mep  (or,  more  generally,  when  M is  nondegenerate)  there  are  2n 
complementary  bases  (counting  multiplicities,  if  necessary).  It  is  also 
well  known  that  binary  n-vectors  can  be  used  as  a way  of  representing  the 
decimal  integers  from  0 to  2n-l  . The  binary  --  i.e.,  base-2  -- 
representation  comes  to  mind  immediately.  One  could,  therefore,  use  the 
set  of  all  binary  n-vectors  appropriately  ordered  as  a way  of  counting 
from  0 to  2n-l  , and  by  the  correspondence  (3)  one  could  use  the 
complementary  bases  in  [I,  -M]  for  this  purpose.  Taking  the  latter 
scheme  and,  in  turn,  its  association  with  complementary  pivot  methods 
into  account,  one  wants  consecutive  integers  to  be  represented  by 
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algorithmically  adjacent  bases.  From  the  standpoint  of  the  methods 
under  consideration  here,  this  means  the  complementary  bases  representing 
consecutive  integers  should  differ  by  exactly  one  column  - or,  equiva- 
lently, that  their  binary  representations  c should  differ  in  exactly 
one  component.  This  being  so,  the  process  of  counting  from  0 to 
2n-l  with  base-2  vectors  cannot  correspond  to  an  execution  of  either 
of  these  complementary  pivot  methods.  For  example,  the  base-2  repre- 
sentations of  2n~^-l  and  2n  ^ differ  in  all  n components. 

The  feature  of  adjacent  complementary  bases  we  wish  to  capture 
with  binary  n-vectors  is  called  the  "unit  distance  property",  a term 
used  in  coding  theory  and  switching  theory.  See  [1].  One  code  having 
the  unit  distance  property  is  known  as  Gray  code  - after  F.  Gray  [11]  - 
also  known  as  the  reflected  binary  code. 

Gray  code  can  be  described  in  various  ways.  For  example,  given 
the  integer 

v e {0,  1,  ....  2n-l } 

one  can  first  write  down  its  base-2  representation* 

bn(v)  = (b-j bn)  (4) 

where 

v = " b.2n_j  . (5) 

j=l  J 


♦At'this  point,  our  notation  Is  somewhat  unconventional,  but  it  serves 
our  purpose  well . 
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Then  the  Gray  code  representation  for  v is  the  binary  n-vector 


9 (v)  ~ ( 9i » •••*  9n) 


whose  components  satisfy  the  relations 


j 91  " b]  (6) 

( 9i  = bi  + bi-i  lmd  2)  , i = 2,  ....  n . 

Table  1 below  contrasts  the  base-2  and  Gray  codes  for  n = 3 . ( 

last  column  of  the  Table  will  be  discussed  a little  later.) 


V 

Base-2  Code  b (v) 

Gray  Code  g (v) 

k3(v) 

0 

(0,  0,  0) 

(0,  0,  0) 

3 

1 

(0,  0,  1) 

(0,  0,  1) 

2 

2 

(0,  1,  0) 

(0,  1,  1) 

3 

3 

(0,  1,  1) 

(0,  1,  0) 

1 

4 

(1,  0,  0) 

(1,1,0) 

3 

5 

(1,  0,  1) 

(1,  1,  1) 

2 

6 

(1,1,0) 

(1,  0,  1) 

3 

7 

(1,  1,  1) 

(1,  0,  0) 

- 

Table  1. 


The  Gray  code  representation  of  the  integers  0,  1,  ....  7 as  binary 
3-vectors  illustrates  two  general  properties: 


Another  feature  of  the  sequence  g (v)  illustrated  in  Table  1 
is  that  every  second  3-vector  is  obtained  from  its  predecessor  by  chang- 
ing the  third  (i.e.,  last)  component.  Again,  this  feature  is  true  in 
general. 

Our  reason  for  devoting  so  much  attention  to  the  Gray  code  is  to 
provide  some  background  for  the  following  assertion  which  is  in  fact  the 
central  point  of  the  paper:  The  complementary  bases  which  arise  in  the 
numerical  examples  of  Murty  and  Fathi  correspond  (via  (3))  to  counting 


in  Gray  code  from  0 to  2-1  . 

After  proving  this  assertion,  we  consider  in  Section  5 the  general 
question  of  what  conditions  on  the  data  must  be  fulfilled  to  make  the 
two  complementary  pivot  methods  follow  a specified  pivot  sequence.  Once 
this  is  done,  the  results  of  Section  4 prove  the  existence  of  data  for 
which  these  conditions  are  satisfied  relative  to  the  specific  pivot  se- 
quence spelled  out  there. 

To  facilitate  the  work  of  Section  4,  we  explore  in  Section  3 some 


properties  of  a function  kn:  {0,  1,  ...,  2 -2}  -Ml,  2 n} 


Definition.  For  all  decimal  integers  v satisfying  0 <_  v <_  2-2  , 


let  kn(v)  denote  the  index  of  the  one  component  in  which  gn(v)  and 


gn(v+l)  differ. 


Thus,  using  the  abbreviation  k = k (v)  we  have 

n 


(7) 


3k(v)  + 9k^v  + ^ = 1 

g^(v)  + g^v  + 1)  = 0 (mod  2),  j / k . 

In  terms  of  complementary  bases,  the  knowledge  of  kn  tells  us  for 
each  v = 0,  1,  ....  2n-2  which  column  of  a given  basis  Is  the  next  to 
be  changed. 

It  is  very  easy  to  verify  that  for  n = 2 we  have 


V 

0 

1 

2 

l<2(v) 

2 

1 

2 

Table  2. 


The  function  k3  has  already  been  displayed  in  Table  1.  Notice  that  the 
corresponding  sequences 


3.  Properties  of  the  function  kp  . 


A comparison  of  k3  with  k^  shows  an  interesting  recursive 


property  of  kn  . 


Proposition  1.  For  all  n > 3 and  0 < v < 2-2 


k (v)  = 


n i f v = 2y 

kn_-j  (y)  if  v = 2p  + 1 . 


Proof.  Let  bn(v)  = (b"(v),  ....  b"(v))  and  gn(v)  = (g"(v),  ....  g”( v) ) 
be  the  base-2  and  Gray  code  representations  of  v , respectively.  There 
are  obviously  two  cases  to  be  considered. 

Case  1 . v = 2y  (i.e.,  v is  even).  This  case  is  trivial  since 
b"(v)  = 0 if  and  only  if  v is  even.  It  follows  that  b"(v  + 1)  = 1 
and  b^v)  = b^v  + 1)  for  j = 1,  ....  n-1  . Hence  by  (6), 

9J(v)  = g"(y  +1)  for  j = 1,  ....  n-1  and  g"(v)  f g"(v  + 1)  . 

This  means  kn(v)  = n . 

Case  2.  v = 2y  + 1 (i.e.,  v is  odd).  We  have  b"(v)  = 1 and 

b^J ( v + 1)  = 0 . Moreover, 


bjV,(v  + 1)  = bJJ_1  (v)  + 1 (mod  2)  . 


By  definition. 


9„(v)  = b"(v)  + b"_j(v)  = 1 + b^_-,(v)  (mod  2) 
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g"(v  + 1)  = b”(v  + 1)  + b”_i(v  + 1)  = 0 + b"  .j(v)  + 1 


(mod  2)  . 


Hence  g"(v)  = gJJ(v  + 1)  whtch  Implies  kn(v)  < n . 
Since 


v = 2p  + 1 = z b?(v)2n-^  and  bn(v)  = 1 
j=l  J n 


we  have 


whence 


which  means 


V - 1 = 2p  = Z b5(v)2n‘j 
j=l  J 


n-1  . . 

v - E b?(v)2n_1"J 
j=l  J 


bn"](p)  = (bJ'V) bj-](p))  = (bj(v),  ....  bJJ_1  (v) ) 


On  the  other  hand 


v + 1 = 2(w  + 1)  = z b"(v  + l)2n"j 

J-l  J 


+ 1 - l b"(v  + l)2n"1-j 
j=l  J 


' :v 


»’■! 


which  means 


bn_1(p  + 1)  = (b"-\v  + 1),  ....  !>;;](„  + 1)) 

= (b"(v  + 1),  ....  bj^_-|(v  + 1))  . 

Now 

gn(v)  = (b"(v),  b2(v)  + blj’(v),  ....  bJJ(v)  + b”_^(v))  (mod  2) 

gn(v+l)  = (b!j*(v+l ) , b^v+l)  + b"(v+l),  ....  b"(v+l)  + bJJ_^(v+l)) 

(mod  2) 

whereas 

gn-1(p)  = (b"(v),  b^(v)  + b"(v),  ....  b^_-|(v)  + bj}_2(v))  (mod  2) 

g11  ^(p+1)  = (b"(v+l),  b2(v+l)  + b"(v+l),  ....  b^_.j(v+l)  + b”_2(v+l)) 

(mod  2) 


It  is  now  clear  that  kn(v)  = k ^ (p)  . □ 

There  is  no  problem  about  evaluating  kn(\>)  when  v is  even. 
For  odd  values  of  v we  can  describe  the  kn(v)  as  follows. 
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Proposition  2.  If  n >_  2 and  0<v<2n-2  Is  an  odd  number,  then 

kn(v)  = n - t 

where  t Is  smallest  positive  Integer  such  that 
v - (2*  - 1)  = 2%^  , Pj.  even. 

Proof.  The  assertion  Is  trivial  when  n = 2 and  v = 1 . Suppose 

n ^ 3 , and  write  v « 2p.  + 1 . If  p^  is  even,  the  assertion  is  true 

with  t = 1 by  Proposition  1.  If  p^  Is  odd  then  p^  = 2p2  + 1 and 
2 2 

v = 2 p2  + (2  - 1)  . If  p2  is  even,  we  have 


kn(v)  = kn-l^l)  = kn-2(u2*  = n ' 2 • 

So  t = 2 . If  p2  Is  odd,  the  process  can  be  repeated  until  we  reach 
an  expression  of  the  form 

v = 2tpt  + (2*  - 1) 

for  which  pt  is  even  and  minimal.  In  this  case  we  have  by  repeated 
application  of  Proposition  1 

kn^  = kn-l^l^  = kn-2^y2^ kn-tW*  = n " 1 ’ 

and  this  is  the  assertion.  □ 


There  Is  a fast  way  to  build  up  kn  from  its  predecessor  kn  1 


if  0 £ v < 2n_1  - 2 (9) 

if  v = 2n_1  - 1 (10) 

+ 1 if  2n_1  < v < 2n  - 2 . (11) 

if  2n_1  < v < 2n  -2  (12) 

1 

and 

kn(2n_1  - 1 - v)  - kn(2n-1  - 1 + v)  if  0 < v < 2n_1  - 1 . (13) 

Proof.  All  these  claims  are  easily  checked  when  n = 3 . Assume,  in- 
ductively,  the  statements  are  true  for  all  m such  that  3 _<  m _<  n - 1 . 

First,  suppose  0 v _<  2n  ^ - 2 . All  such  v belong  to  the  do- 

main  of  kn_^  . If  v is  even,  then  of  course 

kn(v)  = n = (n  - 1)  + 1 = kn-1(v)  + 1 . 

If  v is  odd,  there  exists  an  integer  t as  in  Proposition  2 such  that 


kn(v)  = n - t = (n-l-t)  + 1 = k ,(v)  + 1 . 


Proposition  3.  For  all  n ^ 3 
(kn-l(v)  ♦ 1 

kn(v)  - 

Moreover, 


■■I 

£ 


This  proves  (9). 

Now  let  v = 2n-^  - 1 . Then  v = 2y  + 1 where  w * 2n-2  - 1 . 

It  follows  from  Proposition  2 that 

kn(2n_1  - 1)  = n - (n  - 1)  = 1 

which  proves  (10). 

1 2 

Next,  suppose  2 <_  v <_  2 - 2 . For  all  such  numbers,  we  can 

write 

v = 2n  ^ + v where  0 ^ £ 2n  ^ - 2 . 

This  means  v belongs  to  the  domain  of  kn_^  . If  v is  even,  we  have 

kn(v)  = n . But  In  this  case  “ = v - 2n_1  is  also  even,  whence 

kn_.j(v)  = n - 1 . This  proves  (11)  for  all  even  v . If  v is  odd, 

n 2 n-l 

we  may  write  v = 2y  + 1 where  2n-<:  <_  y <_  2 - 2 . Accordingly,  we 

have 

p = 2n  2 + y where  0^y^2n2-2  . 

By  Proposition  1 

k(v)  = k .(y) 
n n- 1 

k (v)  = k .(77) 
n n- 1 

and  by  Induction, 
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kn-l^  = kn-l^  * 

Putting  these  three  equations  together  we  get  (12).  Using  (12)  and 
induction  on  (11)  we  get 

*„<»>  - *„(■>  - k„.iro  ♦ 1 ■ W*  - 2"_1)  * ' 

which  gives  (11)  for  n . 

Finally,  equation  (13)  asserts  that 

kn(0),  kn(l),  ....  kn(2n  - 3),  kn(2n  - 2)  (14) 

is  a palindrome.  This  has  already  been  seen  for  n = 2 (and  n = 3 ) 
and  since  each  sequence  (14)  is  obtained  from  its  predecessor 

kn-l(zn"’  - 3)'  Vl'2"'’  - 21 

by  putting  the  number  n at  the  beginning  and  end  of  the  sequence  as 
well  as  between  every  consecutive  pair  of  its  terms,  it  follows  (by 
induction)  that  (14)  is  indeed  a palindrome.  □ 
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4.  Application  to  specific  numerical  examples. 


: i 
\ 


I 


[ I 
| » 


In  this  section  we  aim  to  prove  that  solving  the  specific  linear 
complementarity  problems  of  Murty  and  Fathi  by  the  complementary  pivot 
methods  of  Lemke  and  Murty  (called  complementary  pivot  methods  I and  II, 
respectively)  gives  rise  to  a special  pivot  sequence,  namely  the  one 
corresponding  to  the  function  kn  studied  in  Section  3. 

4.1  Murty1 s examples. 

To  exhibit  exponential  computational  requirements  with  complemen- 
tary pivot  method  I,  Murty  defines  q = q(n)  and  M = M(n)  as  follows: 

(0  if  i < j 

1 if  i = j (15) 

2 if  i > j 

Remark  1.  For  all  n 2 , the  matrix  M(n)  is  lower  triangular  and  has 
a positive  diagonal.  It  is  therefore  trivial  to  verify  that  M(n)  e p . 
A significant  fact  about  problems  of  this  sort  is  that  there  Is  a natural 
(greedy)  way  to  solve  (q,  M)  for  any  q . For  i = 1,  ....  n sequen- 
tially define 

A ^ 

z.  = max  {0,  -(q^  + I . (16) 

A A 

The  solution  to  (q,  M)  is  then  q + Mz,  z . For  all  of  the  problems 
(q(n),  M(n))  given  by  (15),  It  Is  Immediate  that  z^  = 2n  and 
^ = 0 1=2,  ....  n . 


q.(n)  = - S 2j 
1 j=n+l-1 
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Remark  2.  According  to  the  definition  (15) 


I > 


- 8 

1 0 0 

' 

II 

co 

cr 

-12 

M(3)  = 

2 1 0 

f. 

-14 

2 2 1 

qn(n)  = -(2n+1  - 2)  = -2(2n  - 1)  . 


This  fact  comes  into  play  later  on. 

Remark  3.  Murty  solves  (q(3),  M(3))  where 


in  complete  detail  by  complementary  pivot  method  I.  After  modifying 
this  algorithm  to  its  equivalent  parametric  (principal  pivoting)  form 
(see  [15],  [7],  [13])  one  can  see  that  the  binary  representations 
(via  (3))  of  the  8=2  bases  visited  are  precisely  the  Gray  code 

representations  of  the  numbers  0,  1 7 . One  can  also  see  that 

the  pivot  sequence  is  k^  , that  is  3,  2,  3,  1,  3,  2,  3 . Analogous 
statements  for  (q(2),  M ( 2 ) ) are  even  easier  to  check,  and  in  this 
case  one  can  draw  a figure  which  nicely  illustrates  the  solution  process. 

Our  goal  now  is  to  generalize  the  observations  made  in  Remark  3 
to  (q(n),  M(n))  for  all  n 3 . Formally,  we  have 


Theorem  1.  The  parametric  form  of  complementary  pivot  method  I applied 
to  (q(n),  M(n))  as  defined  in  (15)  will  generate  the  pivot  sequence 
kn(0),  kn(l),  ....  kn(2n  - 2)  . 


Proof.  The  assertion  is  true  for  n = 2 (as  well  as  n = 3 ).  Assume 


It  is  true  for  n - 1 2 and  consider  the  problem  (q(n),  M(n))  . 

An  additional,  easily-proved  fact  is  that  when  n = 2 , the  cri- 
tical values  of  zQ  are  z°  = 0 , Zq  = 2 , Zg  = 4 , and  zjj  = 6 . 

In  other  words,  Zg  = 2v  . (The  same  can  be  shown  when  n = 3 .) 

We  adopt  this  as  another  part  of  the  inductive  hypothesis. 

For  (q(n) , M(n))  let  zQ  = 2n+1  - 2 = 2(2n  - 1)  . Then  define 

e(n)  = (1 1)T  e Rn 

and 

q(n)  = q(n)  + zQe(n)  . 

It  is  very  easy  to  verify  that 

2n  - 2 

q(n)  = _ . (17) 

q(n-l ) 

The  task  is  now  to  solve  (q(n)-Zge(n) , M(n))  for  all  zQ  € [o,  TQ]  . 
To  do  this,  we  write 

[0,  2n+1-2]  = [0,  2n-2]  U (2n-2,  2n)  U [2n,  2n+1-2]  . 

Now  obviously  for  all  zQ  G [0,  2n-2]  we  must  have  z1  = 0 in 
the  solution  of  (q(n)-Zge(n) , M(n))  . Looking  at  the  structure  of 

1 0 

M(n)  = 

2e(n-l ) M(n-l) 
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and  taking  note  of  (17),  we  see  that  for  Zg  G [0,  2n-2]  we  must  be 
solving  (q(n-l),  M(n-l)  by  the  same  method.  By  the  inductive  hypo- 
thesis we  obtain  the  pivot  sequence  kfi_1(v)  for  v = 0,  1,  ....  2n_1  - 2 . 
For  (q(n),  M(n))  , this  (partial)  pivot  sequence  translates  by  (9)  into 

kn(v)  = kn-l(v)  + 1 for  v = 0,  1,  ....  2n  1 - 2 . 

Next  we  note  that  for  all  Zg  G [2n”^~2,  2n-2]  the  solution  of 
(q^(n)-Zge(n) , M(n))  must  have  z^  > 0 and  w1  = 0 . (See  Remark  1.) 

By  making  z^  basic  in  place  of  w1  , we  have  a pivot  sequence  which 
agrees  with  kn  up  to  v = 2n-^  - 1 . For,  by  (10)  kn ( 2n“ 1 -1 ) = 1 . 

Thus, 

gn(2n’l)  = (1,  1,  0,  ....  0) 

2n-l 

corresponds  to  the  complementary  basis  B . This  is  the  complemen- 
tary basis  for  (q(n)-Zge(n) , M(n))  for  all  zQ  G (2n_1-2,  2n)  . 

We  now  want  to  solve  the  problem  for  zQ  G [2n,  2n+^-2]  . Since 
z^  > 0 and  w^  = q^n)  - zQ  + z^  =0  when  zQ  is  in  this  interval, 
we  have 

Z]  = z0  - (?n  - 2)  (18) 

and  substituting  for  z^  via  (18),  the  values  of  w2»  ...»  wn,  z2, 

....  zn  must  solve 

(q(n-l)  + 2(z0-2n+2)e(n-l)  - Zge(n-l),  M(n-l)) 


-20- 


(19) 


r 


| ! 


i 


for  all  Zg  e [2n,  2n+^-2]  . This  can  be  written  as 

(q(n-l)  + 2e(n-l ) - (2n+1 -2-zQ)e(n-l ) , M(n-l))  . 

Notice  that  for  Zg  £ [2n,  2n+^-2]  we  have 

2n  - 2 >_  (2n+1  -2-zQ)  >0  . 

So  (19)  is  just 

(q(n-l)  + 2e(n-l ) - se(n-l),  M(n-l))  , 2n  - 2 > c > 0 . (20) 


In  solving  (20),  we  get  the  same  sequence  of  complementary  bases  as 
in  solving  (q"(n-1)  - zQe(n-l),  M(n-l))  for  zQ  e [0,  2n-2]  except 
that  they  occur  in  reverse  order.  By  (13)  and  the  inductive  hypothesis 
we  obtain  the  remainder  of  the  sequence  kn  . □ 

This  argument  proves  Murty's  result:  There  are  2n-l  basis  changes 
when  the  parametric  form  of  complementary  pivot  method  I is  applied  to 
(q(n),  M(n))  . It  also  shows  that  the  complementary  bases  encountered 
correspond  to  the  Gray  code  representations  of  the  numbers  v = 0,  1 , 
....  2n-l  . As  noted  by  Murty,  the  critical  values  of  the  parameter 
Zq  are  of  the  form  2v  for  v = 0,  1,  ....  2n-l  and  the  half  line 
of  points  q(n)  - Zge(n)  where  0 <_  Zg  <_  2n+^-2  passes  through  all 
2n  complementary  cones. 


-21- 


Ei 


j 

Is 


Murty  applies  complementary  pivot  method  II  to  the  problem 
(-e(n) , M(n) ) . 

Theorem  2.  Complementary  pivot  method  II  applied  to  (-e(n),  H(n)) 
with  M(n)  as  defined  in  (15)  will  generate  the  pivot  sequence 
kn(0),  kn(l),  ....  kn(2n-2)  . 

Proof.  This  is  easily  verified  in  the  case  where  n = 2 . Now  suppose 
it  is  true  for  n - 1 > 2 . We  have 


' -1  ' 

1 0 

e(n)  = 

and  M(n)  = 

-e(n-l) 

2e(n-l)  M(n-l) 

By  the  way  the  method  is  defined,  we  must  have  = 0 (non  basic) 
until  (-e(n-l),  M(n-l))  is  solved.  By  the  inductive  hypothesis, 
this  generates  the  pivot  sequence  kn_^v)  for  \>  = 0,  1,  ...,  2n-1-2  . 
At  the  termination  of  the  latter  solution  process  we  have  - in  terms 
of  the  problem  (-e(n),  M(n))  - the  variable  z^  basic  and  all  other 
zi  non  basic.  The  next  variable  to  become  basic  is  z^  , so  by  Propo- 
sition 3 we  have  the  pivot  sequence  kn(v)  for  0 < v < 2n-l  . Now 
with  z.j  and  z^  basic  and  z.  non  basic  for  3 <_  i <_  n we  must 
once  again  have  to  solve  the  problem  (-e(n-l),  M(n-l))  . (We  remark 
that  w.j  will  never  become  basic  again.)  By  the  inductive  hypothesis, 
this  requires  the  sequence  k^  ^(v)  for  0 <_  v £ 2n"^-2  . Since  this 
sequence  is  a palindrome,  it  follows  that  the  entire  pivot  sequence 
for  solving  (-e(n),  M(n))  is  just  kn  . □ 


mentary  bases  used  correspond  to  the  Gray-code  representations  of  the 
numbers  0,  1 , ....  2n-l  . 

4.2  Fathi's  examples. 

The  main  point  of  Fathi's  paper  is  that  there  are  symmetric  P- 
matrices  (which  of  course,  must  therefore  be  positive  definite)  for 
which  the  same  computational  effort  is  required  as  in  the  solution  of 
Murty's  examples  by  the  same  methods.  Indeed,  Fathi  constructs  a 
special  matrix  out  of  Murty's  matrix  M(n)  . Before  defining  this 
matrix,  we  wish  to  call  the  reader's  attention  to  the  fact  that  there 
are  some  significant  notational  differences  between  the  following 
development  and  what  is  to  be  found  in  Fathi's  paper.  So,  to  begin 
with,  let 


F(n) : = M(n)M(n) 


(21) 


This  matrix  is  obviously  symmetric  and  positive  definite.  Furthermore, 
it  has  the  following  suggestive  structure: 


F(n)  = 


1 2e(n-l ) 

2e(n-l ) F(n-l)  + 4E(n-l)J 


(22) 


where 


E(n) : = e(n)e(n) 
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Fathi's  analysis  (and  ours)  makes  use  of  two  related  rmtrices,  namely 


F(n) : = F(n)  + 4E(n) 


(23) 


and 


F(n):  = F(n)  - 4G(n) 


(24) 


where 


G(n) : = [0,  e(n)e(n-l )T]  . (25) 

We  notice  at  once  that  F(n-l)  is  a principal  submatrix  of  F(n)  . 

For  reasons  of  simplicity,  we  reverse  the  order  of  the  complementary 
pivot  methods  considered  and  look  at  the  solution  of  Fathi's  problem 
(-e(n),  F(n))  by  means  of  complementary  pivot  method  II. 

Examining  Fathi's  solution  of  (-e(3),  F(3) ) we  see  that  comple- 
mentary pivot  method  II  generates  the  pivot  sequence  k^  . One  can 
easily  verify  the  analogous  statement  (-e(2),  F(2))  . In  Theorem  3 
we  prove  this  in  general.  The  proof  (naturally an  inductive  one)  makes 
critical  use  of  an  observation  of  Fathi's.  We  state  this  without  proof 
as 

Lemma  1.  For  all  n ^ 2 , the  pivot  sequences  generated  by  complementary 
pivot  method  II  in  the  solution  of  (-e(n),  F(n) ) , (-e(n),  F(n))  , 

and  (-e(n),  F(n))  are  the  same. 

Now  we  have 
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Theorem  3.  Complementary  pivot  method  II  applied  to  (-e(n),  F( n) ) 
will  generate  the  pivot  sequence  kn(0),  ....  kn(2n-2)  . 

Proof.  The  statement  is  true  for  n = 2 , so  assume  it  is  true  for 
n - 1 2 . To  solve  (-e(n),  F(n))  one  must  first  solve  the  sub- 

problem (-e(n-l),  F(n-l))  before  can  be  made  basic.  By  the 
inductive  hypothesis  and  Lemma  1 we  see  that  the  solution  of 
l-e(n),  F(n) ) by  complementary  pivot  method  II  begins  with  2n~^-l 
pivots  given  by  kn(v)  for  0 <_  v £ Zn~  -2  . The  next  pivot  must 
be  the  ex:hange  of  w^  and  . Hence  we  have  agreement  with 
kn(2n-1-l)  . When  this  stage  is  reached,  z^ , z^,  w^,  ....  wn 
are  basic  variables.  The  corresponding  tableau  has  the  form 


1 

W1  w2  z3  zn 

Z1 

3 

5-22  •••  2 

Z2 

-1 

-2  1 -2  ...  -2 

W3 

-1 

CM 

CM 

1 

• 

• 

• 

• 

F(n-2) 

• • 

• 

w , 

• 

-1 

• • 

-2  2 

w 

According  to  the  ruler  of  the  method,  we  must  solve  the  subproblem  cor- 
responding to  the  indices  2,  3 n . This  is  easily  recognized 

as  (-e(n-l),  F(n-l))  . By  the  inductive  hypothesis  and  Lemma  1,  we  see 
that  pivot  sequence  for  this  part  of  the  problem  is  just  the  last 
2n"^-l  terms  of  the  sequence  kR  , i.e.,  kn(v)  for  2n~^  < v < 2n-2  . 
After  these  pivots,  the  problem  is  solved.  Hence  complementary  pivot 
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method  II  applied  to  (-e(n),  F(n))  generates  the  pivot  sequence 
kn(v)  for  0 < V < 2n-2  . □ 

To  obtain  a computational  complexity  result  for  complementary 
pivot  method  I,  Fathi  uses  the  matrix  we  have  called  F(n)  , but  instead 
of  q(n)  as  defined  in  (15),  he  considers  a more  general  class  of 
vectors.  The  central  idea  in  this  instance  is  to  show  the  existence 
of  a line  L(n)  which  passes  through  all  2n  complementary  cones  de- 
fined relative  to  F(n)  . 

The  construction  goes  as  follows.  Let  f>(n)  e Rn  satisfy 

0 < Pn(n)  < — <P1(n)  . (26) 

Note  that  p^(n)  e int  . For  purposes  of  the  proof,  it  would  suffice 
to  let  P^n)  equal  0 in  (26).  The  desired  line  is 

L(n):  = (p(n)  - zQe(n):  z^R}  . (27) 

A figure  corresponding  to  the  parametric  linear  complementarity 
problem  (p(2)-zQe(2) , F(2))  shows  that  L ( 2 ) passes  through  the  four 
complementary  cones  defined  relative  to  F(2)  . It  also  shows  that  if 
Zg  is  suitably  large  and  positive,  then  p (2 ) : = ”p(2)  - Zge(2)  lies 
in  the  interior  of  pos  [-F  .J , * Moreover,  the  line  segment  be- 

tween p(2)  and  p\2)  meets  all  four  complementary  cones  relative  to 
F(2)  . Indeed,  solving  the  linear  complementarity  problem  (p(2),  F( 2 ) ) 
by  complementary  pivot  method  I in  parametric  form  yields  the  pivot 
sequence  k 2 . 
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As  in  the  previously  considered  cases,  this  fact  points  to  an  in- 
ductive proof  of  the  natural  generalization  of  this  observation.  We 
shall  not  give  a formal  proof  of  the  following  result,  but  simply 
remark  that  using  the  obviously  true  case  of  n = 2 as  an  inductive 
hypothesis  with  Fathi's  own  inductive  proof  of  his  Theorem  3.4  we  have 

Theorem  4.  If  p(n)  € L(n } lies  in  the  interior  of  the  complementary 
cone  pos  [-F^,  I 2 I.n]  , then  the  parametric  version  of  comple- 

mentary pivot  method  I applied  to  (p(n),  F(n) ) will  generate  the 
pivot  sequence  k^fO),  ....  kn(2n-2)  . 

{ 

Thus,  the  two  complementary  pivot  methods  applied  to  the  linear 
complementarity  problems  set  down  by  Murty  and  Fathi  give  rise  to  comple- 
mentary bases  which  correspond  to  the  Gray-code  representations  of  the 
integers  0,  1,  ....  2n-l  . 


I 
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5.  Conditions  for  achieving  a prescribed  pivot  sequence. 

In  this  section,  we  are  concerned  with  the  relationship  between 
the  data  q and  M of  a linear  complementarity  problem  and  the  pivot 
sequence  followed  in  the  solution  of  (q,  M)  by  complementary  pivot 
methods  I and  II. 

Throughout  this  section  we  assume  that  M e Rnxn  js  a P-matrix. 

We  know  then  that  for  all  q £ Rn  , the  problem  (q,  M)  has  a unique 
solution  [17]  and  either  of  the  aforementioned  methods  will  find  it. 

The  "results"  of  this  section  serve  two  purposes.  First,  they 
provide  a brief  review  of  two  complementary  pivot  methods  under  con- 
sideration here.  (This  may  be  helpful  for  those  not  already  conversant 
with  them.)  Second,  and  no  less  important,  they  suggest  the  possibil ity 
of  constructing  other  nasty  problems,  about  which  more  will  be  said  in 
Section  6. 

Both  of  these  complementary  pivot  methods  use  an  algebraic  opera- 
tion know  as  principal  pivoting.  In  general  (i.e.,  for  M not  neces- 
sarily in  P ),  if  m....  f 0 , the  principal  pivotal  transform  of  M ob- 
tained by  pivoting  on  m...  is  the  matrix  M'  with  elements  m'.j 
satisfying  the  rules 

mii  ■ mH  ) 

mij  = "miimi j / j * 1 

mjf  = mjA'  ) 

mj* = mjt  ■ • 3 * 1 * 1 • 

A salient  point  in  the  present  discussion  is  a result  due  to  Tucker  [19] 
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which  states  that  if  M £ P then  so  is  every  principal  pivot;  1 trans- 
form M1  . Thus  the  assumption  that  H G P guarantees  that  all  the 
indicated  pivot  operations  are  possible. 

5.1  Complementary  pivot  method  I (Lemke  [13]). 

In  the  present  circumstances,  Lemke's  method  can  be  put  into  para- 
metric form  as  follows.  Given  a problem  (q,  M)  of  order  n with 
q f_  0 (for  nontriviality)  define  the  positive  real  number 

z0  = - min  q 
llifn  1 

and  let  e = (1 l)^6Rn  , Then 

q":  = q + Ige  _>  0 . 

(Recall  what  was  done  in  the  proof  of  Theorem  1.)  To  solve  (q,  M) 
we  may  consider  the  parametric  linear  complementarity  problem 

(q-zQe,  M)  , 0<z0<z0  . (28) 

For  all  zQ  g [o,  zQ]  , the  problem  (q-Zge,  M)  has  a unique  solution. 
For  Zg  = 0 , the  solution  is  obvious.  What  we  want  is  the  solution  for 
Zg  = Zg  . This  can  be  thought  of  as  the  homotopic  approach  to  the 
problem. 

As  Zg  runs  from  0 to  zQ  , the  vector  c[  - zQe  passes  through 
several  complementary  cones  along  the  line  segment  from  q to  q . 

The  critical  values  Zg  of  Zg  correspond  to  points  where  this  line 
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segment  meets  the  boundaries  of  the  complementary  cones.  For  example, 
Zq  = 0 is  the  largest  value  for  which  q - zQe  e pos  B°  where 
B°  = I . Thus  z^  = 0 is  a critical  value  of  zQ  . For  any  larger 
value  of  Zq  , c[  - zQe  belongs  to  a different  complementary  cone. 

In  the  nondegenerate  case,  this  cone  is  uniquely  determined.  So,  under 
this  assumption,  for  the  critical  value  Zq  , there  will  be  a unique 
index  k(0)  such  that 

‘ z0e^k(0)  = 0 * 

In  the  degenerate  case,  several  components  of  q - z^e  may  vanish. 

One  can  make  the  choice  unique  and  the  method  finite  by  using  a least- 
index  rule  [3],  We  shall  say  more  about  this  later;  for  the  moment  we 
retain  the  nondegeneracy  assumption. 

According  to  (29),  the  line  segment  from  q to  q first  enters 
the  complementary  cone  pos  B^  where 

B1  = [I.-, , ...,  I.k(0)-r  -M.k(0)’  1 -k(0)+l » •••»  !-n]  ‘ 

The  process  then  continues  after  the  system  is  put  into  canonical  form 
with  respect  to  the  complementary  basis  B^  . This  corresponds  to  a 
principal  pivot  on  m^Qj  • Equivalently,  we  take  the  equation 

Iw  - Mz  = q - zQe 

and  multiply  through  by  ( b"* )“"*  to  obtain  the  up-dated  system.  We 
must  next  determine  the  critical  value  Zq  , that  is,  the  largest  value 
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of  zQ  for  which  q - zQe  g pos  B . However,  in  this  and  all  sub- 
sequent stages,  care  must  be  taken  to  make  sure  that  the  next  critical 
value  does  not  exceed  Zg  . In  particular  if  Zg  <_  zQ  < Zg+^  , then 
q G pos  Bv  and  the  process  stops  as  soon  as  the  easily  found  solution 
of  (q-I  e,  M)  = (q,  M)  is  written  down. 

Suppose  the  solution  process  requires  exactly  i pivots  (basis 

changes).  This  means  sequences  of  i + 1 complementary  bases 

Oi  p o 1 a 

B , B , ....  B and  critical  values  zQ,  Zg,  ...,  Zg  are  determined 

in  the  manner  sketched  above.  Each  basis  is  obtained  from  its  prede- 
cessor by  a simple  column  change.  For  v = 0,  1 , ... , fc-1  let  k(v) 
denote  the  index  of  the  column  in  Bv  to  be  exchanged.  Thus 


I 


BV*J  x « (I  \ , -M  \ }\  {Bv.  / 

• k( v)  *k(v)’  *k(v)  *k(v) 


Bvt]  = Bv.  , 
•J  *J  * 


j t k ( v) 


If  we  have  Bu  = I , and  we  know  the  pivot  sequence  (function) 
k(0),  k(l),  ....  k(Ji-l)  , then  we  know  all  the  complementary  bases 
used  in  the  process. 

What  we  seek  is  the  connection  between  the  data  of  (q,  M)  and  a 
given  pivot  sequence  k . The  question  is:  When  will  complementary 
pivot  method  I (in  parametric  form)  generate  a particular  sequence  of 
pivots? 

This  question  has  an  obvious  answer  in  terms  of  some  determinantal 
inequalities.  To  facilitate  their  expression,  we  introduce  the  following 
notation.  Suppose  A is  an  m*n  matrix,  b is  an  m-vector,  and 
k G (l , ....  n}  . Then 


A ( k ; b)  i = fA^ , •••»  k»  A ^ * •••»  A#^] 

is  just  the  matrix  A but  with  b in  place  of  A^  . Note  that  if 
m = n and  A is  nonsingular,  then  by  Cramer's  rule  for  solving  the 
linear  algebraic  system  Ax  = b we  have 

x^  = det  A(k;  b)/det  A , k = 1 , ....  n . 

Concerned  as  we  are  with  complementary  bases  relative  to  [I,  -M] 
where  Hep  , it  should  be  clear  that  given  such  a basis,  B , its 
determinant  is  positive  or  negative  according  to  whether  B contains 
an  even  or  odd  number  of  columns  of  -M  . Moreover,  given  B°  = I 
and  a pivot  sequence  k(0),  ....  k(t-l)  , the  complementary  bases  B 
it  generates  will  contain  an  even  or  odd  number  of  columns  from  -M 
according  to  whether  v is  even  or  odd.  Thus 

(-l)v  det  Bv  > 0 , 0 £ v < z . 

We  use  this  fact  in  the  proof  of  the  theorem  below. 

Theorem  5.  If  M € p , and  the  solution  of  (q,  M)  by  (the  parametric 
form  of)  complementary  pivot  method  I generates  the  pivot  sequence 

k(0) k(t-l)  and  corresponding  complementary  bases  B°,  ...,  B5,  , 

then  for  v = 0,  ....  t-1 
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(i) 


(-l)v  det  Bv(k(v) ; e)  > 0 ; 


(ID 


(-l)v  det  Bv(k(v) ; q)  < 0 ; 


(f-j-j)  det  i = min  (det  B (j;  q) 

det  Bv(k(v);  e)  1_<j<n  (det  Bv(j;  e) 


(-1 )vdet  Bv(j ; e)  > 0 


and  for  v = *. 


( i v ) either 

max  (-l)adet  B^(j;  e)  ^ 0 
1 <j_<n 


or  else  there  exists  an  index  k(«.)  for  which 
(-l)^det  Ba(k(^) ; e)  > 0 , (iii)  holds  for  v = i and 
(-l)^det  B£(kU);  q)  0 . 


k. 


Conversely,  if  (i)  - (i v)  hold  and  for  all  zQ  e [0,  Zg]  there  are  at 
most  n + 1 zero  components  in  the  complementary  basic  solution  to 
(q  - zQe,  M)  , then  the  sequence  of  pivots  must  be  given  by  the  integers 
k(0),  ....  kU-1)  . 

Proof.  Suppose  the  pivot  sequence  k(0),  ....  k(fc-l)  is  generated. 

Then  for  each  v = 0,  t-1  , the  k(v)-th  basic  variable  in  the 

v-th  basic  set  of  variables  must  decrease  as  zQ  increases.  To  simpli- 
fy the  notation,  let  k = k(v)  , and  let  wjj  be  the  k-th  basic  vari- 
able relative  to  the  v-th  basis  Bv  . Then 
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3W 

3Z 


* = ((Bv)-1(-e)).  = det.Bv(.k-,  -e)  , 0 # 


det  B' 


H 


i 


i * 


Since  M G p , we  have  (-1)vdet  Bv  > 0 , the  Inequality  (1)  now  follows. 

Next,  consider  the  calculation  of  the  critical  value  Zg  . The 
up-dated  right-hand  side  of  the  equation 

Iw  - Mz  = q - zQe 

is  of  the  form 


(Bv)_1q  - zn(Bv)"1e  . 


The  critical  value  Zg  makes  the  k-th  component  of  the  vector  in 
(31)  equal  zero.  Hence  for  each  v = 0 t-1 


z„  = 


,V\-1  _\  / /r>V\“l 


( (B  ) *e) 


k(v) 


«B  ) '•>«*) 


«BV>  *)«„,  0 


However,  for  these  values  of  v , Zg  <_  Zg  . The  inequality  ( 1 i ) now 
follows  from  (i).  The  fact  that  Zg  is  the  critical  value  for  the  Bv 
implies  (iii). 

When  the  last  basis  B*  is  obtained,  it  must  be  the  case  that 
either  (B*)“^(-e)  > 0 in  which  case 
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0 j?-l 

and  none  of  the  variables  w^.  decreases  as  Zg  increases  from  Zg 

to  zQ  or  else  (Be)~^(-e)  has  a negative  component  and  it  is  possible 
to  define  a new  critical  value  Zq>.7q  as  above.  In  particular,  there 
must  exist  an  index  k(t)  for  which  (iii)  holds  and 

(-l)ldet  B^kU);  q)  >_  0 . 

Since  (B£)-1(q  - zQe)  >_  0 for  all  zQ  £ [Zq-1 , zQ]  , it  is  not  neces- 
sary to  execute  another  pivot;  hence  k ( «. ) is  not  part  of  the  pivot 
sequence. 

The  proof  of  the  converse  is  clear  form  the  nondegeneracy  assump- 
tion. Indeed  the  pivot  sequence  is  uniquely  determined  by  the  inequali- 
ties (i)  - (iv).  □ 

Remark.  The  necessary  conditions  stated  in  Theorem  5 are  valid  even 
in  the  degenerate  case,  but  in  that  situation,  the  statement  of  the 
converse  must  be  modified.  If  we  use  a least-index  rule  to  choose  the 
pivot  element,  then  (i)  - (iv)  are  sufficient  provided  the  k(v)  are 
the  least  integers  (indices)  for  which  those  conditions  hold. 

5.2  Complementary  pivot  method  II  (Murty  [15],  [16]). 

In  [15]  and  [16],  Murty  develops  variants  of  the  so-called  Bard 
method*  for  solving  (q,  M)  where  M £ P . The  Bard  procedure  is  a 


*This  sort  of  method  is  related  to  an  earlier  one  of  Zoutendijk  [20, 
pp.  83-90]. 


type  of  principal  pivoting  method  which  starts  with  the  equation  (1) 
and  the  complementary  basis  B°  = I . At  the  v-th  step,  one  has  a 
complementary  basis  Bv  and  the  vector 

V /„Ui -1 

q : = (B  ) q . 

If  qv  21  0 , the  solution  has  been  formed.  Otherwise  any  index 
k = k(v)  such  that  q^  c 0 is  chosen,  a new  complementary  basis 
Bv+1  is  determined  according  to  (30),  and  by  principal  pivoting  (or 
up-dating)  one  has  the  new  vector 


v+1  /„v+l \~1 

q : = (B  ) q . 


At  this  point,  the  steps  outlined  above  are  repeated. 

It  is  known  that  in  the  degenerate  case  this  arbitrary  way  of 
selecting  k(v)  can  lead  to  circling  (alias  cycling):  the  repetition 
of  a sequence  of  complementary  bases.  The  essence  of  Murty's  variant 
is  a specific  rule  for  choosing  k = k(v)  which  makes  the  method 
finite.  In  [14],  he  uses  the  rule 


k(v)  = arg  min  { k : q^  < 0} 


whereas  in  [16]  he  uses  the  largest-index  rule 


k(v)  = arg  max  { k : q^  < 0}  . 
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Under  (32)  or  (33)  or  indeed  under  any  other  specific  rule,  k becomes 
a function  of  v and  the  data  q and  M . 

So,  complementary  pivot  method  II  is  Bard's  method  with  the  largest 
index  pivot  selection  rule  (33).  It  seems  that  Murty's  motivation  for 
using  (33)  rather  than  (32)  in  [16]  was  to  exhibit  the  required  exponen- 
tial computational  effort  of  both  complementary  pivot  methods  with  the 
same  matrix.  Fathi  [7]  also  uses  (33)  as  part  of  complementary  pivot 
method  II. 

Theorem  6.  If  M e P n Rn*n  , the  solution  of  (q,  M)  by  complementary 

pivot  method  II  generates  the  pivot  sequence  k(0),  k(l),  ....  k()i-l) 

0 1 a 

and  corresponding  complementary  bases  B , B , ....  B if  and  only  if 

(i)  for  all  v = 0,  1 , ...»  t-1 

min  C-l ) vdet  Bv(k;  q)  < 0 

l-iKi11 

and 

k(v)  = arg  max  {k:  (-l)vdet  Bv(k;  q)  < 0}  ; 

(ii)  for  v = z 

(-l)*det  B*(k;  q)  > Q , k = 1 n . 

Moreover,  the  bound  i <_  2n  - 1 is  sharp. 

Proof.  Given  a complementary  basis  Bv  we  have 

= ((Bv)-1q)k  = det  Bv(k;  q)/det  Bv  . 
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Since  (-l)vdet  Bv  > 0 for  all  complementary  bases  generated  by  comple- 
mentary pivot  method  II,  conditions  (i)  and  (ii)  are  just  a paraphrase 
of  the  tests  of  the  algorithm.  The  bound  z <_  2n  - 1 is  a consequence 
of  the  fact  that  there  are  2n  complementary  bases  which  (by  Murty's 
finiteness  proof)  cannot  be  repeated  in  the  solution  process.  The 
sharpness  of  the  bound  is  shown  by  Murty's  problem  (-e(n),  M(n) ) . □ 


I 


6.  Some  connections  with  other  literature. 


In  section  4,  we  established  the  correspondence  between  the  numerical 
examples  of  Murty  and  Fathi  and  the  Gray  code  representations,  gn(v)  , of 

the  integers  v = 0,  1 2n-l  . This  sequence  describes  a hamil- 

tonian  path  on  the  unit  n-code. 

As  may  be  expected,  there  is  a sizable  literature  on  Gray  codes 
and  related  matters.  We  mention  just  a little  of  it  here.  Some  papers 
are  concerned  with  Gray  codes  in  bases  other  than  2 and  conversions 
between  number  systems.  See  [4],  [6],  and  [8],  Other  papers  treat  the 
problem  of  describing  paths  on  n-cubes.  See  [10],  for  example.  The 
latter  are  related  to  "snake-in-a-box"  problems  and  codes.  The  Intro- 
duction (in  Volume  I)  and  Bibliography  (in  Volume  II)  of  the  book  [14] 
by  MacWilliams  and  Sloane  is  a rich  source  of  references  on  this  litera- 
ture. See  also  [2], 

By  far  the  most  amusing  article  on  the  subject  of  Gray  codes  is 
the  one  by  Martin  Gardner  [9]  who  points  out  that  the  binary  Gray  code 
is  the  key  to  the  solution  of  the  Chinese  ring  puzzle  and  the  well- 
known  Tower  of  Hanoi  problem.  To  these  may  be  added  a recently-mar- 
keted puzzle  called  "The  Brain."  It  requires  the  solver  to  move  a set 
of  eight  rod?;  from  one  extremal  position  (e.g.,  all  "out")  to  the  other 
(all  "in").  The  puzzle  is  constructed  in  such  a way  that  at  any  stage, 
only  two  of  the  rods  will  move.  It  takes  170  moves  to  solve  the  puzzle, 
for  (in  the  notation  of  Section  2) 

g8(o)  = (o,  o,  o,  o,  o,  o,  o,  o) 


and 
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g8 ( 1 70)  = (1,  1,  1,  1,  1,  1,  1,  1)  . 

Remembering  that  the  papers  by  Murty  and  Fathi  are  concerned  with 
the  computational  complexity  of  complementary  pivot  methods  and  that 
problems  requiring  2n  - 1 pivots  are  exhibited,  one  may  ask  whether 
these  have  a connection  with  the  Klee-Minty  result  [12]  on  the  simplex 
method  of  linear  programming.  In  their  paper,  Klee  and  Minty  devise 
a class  of  linear  programs  in  which  the  feasible  region  is  a specially 
perturbed  unit  n-cube  . The  objective  function  is  just  one  coordinate, 
and  the  simplex  method  with  the  customary  pivot  selection  rules  runs 
through  2n  - 1 pivots  in  finding  the  solution.  In  so  doing,  it 
generates  a hamiltonian  path  cn  the  "cube"  which  can  be  associated  in 
a simple  way  with  the  path  that  arises  from  the  Gray  code.  Their  exam- 
ple [12,  p.  163]  with  e = 0 and  the  coordinates  of  the  vertices  written 
in  reverse  order  illustrates  this  point.  It  is  not  clear  (to  the  author 
at  least)  that  a stronger  link  between  these  phenomena  can  be  found. 

For  the  combinatorial  theorist,  an  appealing  question  is:  What 
"solution  paths"  on  the  unit  n-cube  can  be  realized  - through  the 
correspondence  (3)  - when  a linear  complementarity  problem  (q,  M) 
with  H ep  is  solved  by  methods  of  the  type  discussed  here? 

Very  recently,  Stickney  and  Watson  [18]  have  published  a paper 
focusing  attention  on  Bard-type  algorithms  for  the  P-matrix  case.  For 
a given  problem  (q,  M)  they  define  a directed  graph  G(M/q)  whose 
vertices  correspond  to  c pi ementary  bases  of  [I,  -M]  . The  graph  has 
a directed  edge  from  vertex  A to  vertex  B if  there  exist  an  index 
k € {1 , ...»  n}  such  that 
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(1) 

(A_1q)k  < 0 

(34) 

(ii) 

A-k  G{I-k*  -M.k}VB-k 

(Hi) 

A.j  - B.j  . j t k . 

Clearly,  the  vertices  of  G(M/q)  correspond  to  those  of  the  unit  n-cube 
via  (3).  The  edges  of  G(M/q)  being  defined  by  the  conditions  (34)  are 
just  those  of  the  unit  n-cube  apart  from  their  orientation.  That  is, 
when  M £P  , the  total  degree  of  a vertex  In  G(M/q)  is  n . See 
[18,  p.  324]. 

One  of  the  main  results  of  the  Stickney-Watson  paper  is  that  when 
H ep  , a linear  complementarity  problem  (q,  M)  with  q belonging  to 
a complementary  cone  spanned  by  k columns  of  -M  and  n - k columns 
of  I can  be  solved  by  a sequence  of  k Bard-type  pivots.  Unfortunately, 
in  the  general  case,  one  still  has  no  clue  as  to  which  Bard-type  pivots 
ought  to  be  executed.  Perhaps  future  studies  will  shed  more  light  on 
this.  Clearly,  for  this  insight  to  have  practical  algorithmic  value, 
it  will  be  necessary  to  identify  the  "right"  pivots  without  an  excessive 
amount  of  auxiliary  work. 

Stickney  and  Watson  also  comment  on  the  problem  of  circling  in 
Bard- type  algorithms  and  about  the  effect  of  scaling  in  this  regard. 

In  view  of  the  relationship  between  complementary  bases  and  vertices 
of  the  unit  n-cube,  and  hence  between  principal  pivoting  algorithms 
and  paths  on  the  unit  n-cube,  it  is  tempting  to  speculate  about  the 
Interpretation  of  these  vertices  (binary  n-vectors)  as  encoded  numbers. 

In  particular,  can  the  sequence  {Bv}  of  complementary  bases  be  chosen 
in  such  a way  that  the  corresponding  sequence  (c(Bv)}  of  binary 


j 


L 
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n-vectors  is  strictly  increasing  in  the  sense  of  some  binary  number 
system?  If  so,  this  might  be  a way  to  avoid  circling. 
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